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1.  Introduction 


Dynamic  fracture  mechanics  is  a  specialized  branch  of  engineering  that  deals  with  the  fracture 
phenomena  for  which  the  role  of  material  inertia  cannot  be  ignored.  Dynamic  fracture  is  a 
subfield  of  fracture  mechanics.  In  general,  the  effects  of  inertia  cannot  be  ignored  when  loads 
are  applied  at  relatively  high  rates  to  a  solid  body  with  an  existing  crack  or  from  a  rapidly 
growing  crack.  Normally,  when  such  bodies  are  subjected  to  surface  tractions,  the  crack  may  or 
may  not  propagate,  depending  on  the  magnitude  and  rate  at  which  the  traction  is  applied.  For 
low  rates,  classical  fracture  mechanics  may  be  used  to  solve  the  problem  in  which  fracture  can  be 
predicted  once  the  stress  level  reaches  a  critical  value.  However,  for  high  rates,  the  crack  may 
propagate  instantaneously,  depending  on  the  magnitude  of  the  stress  wave  and  the  rate  at  which 
it  is  applied.  Therefore,  the  stress  field  around  the  crack  tip  is  not  in  equilibrium,  and  the  inertial 
effects  cannot  be  ignored. 

Unlike  classical  fracture  mechanics  which  is  somewhat  well  established,  dynamic  fracture 
mechanics  is  relatively  new  when  compared  to  other  branches  of  engineering.  Although  some 
advances  have  been  made  over  the  past  few  decades,  the  laws  governing  the  behavior  of  cracks 
when  subjected  to  stress  waves  are  still  in  their  infancy;  therefore,  the  approach  to  analyzing 
such  problems  is  based  on  qualitative  reasoning  and  experience  with  known  solutions  to  specific 
problems  (1 ). 

A  popular  concept  in  the  treatment  of  fracture  initiation  and  propagation  due  to  dynamic  loading 
is  the  stress  intensity  concept.  It  provides  a  basis  for  quantifying  the  resistance  of  materials  to 
the  onset  of  growth  of  an  existing  crack.  The  basis  for  this  approach  can  be  traced  back  to  the 
Irwin  criterion  in  classical  fracture  mechanics.  This  criterion  states  that  the  onset  of  crack 
growth  occurs  when  the  stress  intensity  factor  of  a  material  has  risen  to  a  critical  value  called  the 
fracture  toughness,  Kjc.  The  experimental  detennination  of  Abe  can  be  quite  difficult.  Similarly, 
dynamic  fracture  experiments  are  extremely  difficult  and  expensive  to  execute  because  the 
equipments  involved  require  a  steep  learning  curve  and  high  cost.  For  this  reason,  fewer  than 
100  facilities  around  the  world  pursue  dynamic  fracture  experiments. 

During  the  19th  century,  limited  forms  of  high  strain  rate  experimentations  were  performed  on 
military  ordnances  and  industrial  machines.  However,  no  major  breakthrough  can  be  identified 
as  the  launching  pad  for  dynamic  fracture  mechanics  as  an  area  of  research  until  the  early  1970s 
when  the  importance  of  understanding  crack  propagation  and  arrest  was  recognized  in  both 
engineering  and  earth  sciences.  Some  of  the  ideas  developed  for  fracture  under  equilibrium 
conditions  are  also  useful  for  dynamic  fracture.  The  work  by  A.  A.  Griffith  (2)  is  considered  to 
be  one  of  the  earliest  in  equilibrium  fracture  mechanics  that  is  useful  as  a  quantitative  science  in 
material  behavior.  Plane  elasticity  solution  for  a  crack  of  finite  length  in  a  body  subjected  to  a 
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uniform  remote  tension  in  a  direction  normal  to  the  crack  plane  was  the  basis  of  Griffith’s 
original  work.  Griffith,  when  considering  an  ideally  brittle  elastic  body  containing  a  crack, 
recognized  that  the  macroscopic  potential  energy  of  the  system,  which  consists  of  the  internal 
stored  elastic  energy  and  the  external  potential  energy  of  the  applied  loads,  changes  as  the  crack 
grows.  Furthermore,  Griffith  realized  that  due  to  the  growth  of  the  crack  or  crack  extension,  new 
crack  surfaces  are  created.  Therefore,  he  postulated  that  a  certain  amount  of  work  per  unit  area 
of  crack  surface  must  be  expended  at  a  microscopic  level  to  create  that  surface.  This  work  per 
unit  area  is  usually  characterized  by  assuming  that  a  particular  force-displacement  relationship 
controls  the  reversible  interaction  of  atom  or  molecules  across  the  fracture  plane.  Griffith’s 
energy  condition  is  of  particular  interest  to  engineers  and  scientists  because  it  renders  the  need  to 
examine  in  detail  the  actual  fracture  process  at  the  crack  tip  unnecessary. 


There  were  several  important  works  relating  to  fracture  worth  mentioning  after  that  of  A.  A. 
Griffith.  Sir  Neville  Mott  (3)  proposed  a  theoretical  framework  that  includes  inertial  effects 
during  the  rapid  crack  growth  phase.  At  high  crack  speeds,  Mott  realized  that  inertial  resistance 
of  the  material  to  crack  opening  could  become  significant.  Mott  assumed  that  the  crack  growth 
process  was  time  independent  in  order  to  estimate  the  crack  speed  for  a  specific  loading  system. 
However,  it  was  later  shown  that  Mott’s  assumptions  of  the  dynamic  Griffith  crack  were  not 
valid,  and  therefore  the  resulting  conclusions  were  also  not  valid.  George  Irwin  (4),  concerned 
with  cleavage  fracture  of  structural  steels,  developed  another  important  theoretical  idea  for 
applying  work  and  energy  methods  to  fracture  dynamics.  Although  Irwin  adopted  Mott’s 
postulate  of  total  energy,  which  includes  a  work  of  fracture,  he  proposed  that  the  work  of  fracture 
may  be  approximately  represented  as  the  sum  of  two  proportional  terms.  These  are  the  area  of 
fracture  surface  and  the  volume  of  material  affected  by  plastic  deformation.  This  idea  therefore 
expanded  the  applicable  range  of  the  energy  balance  fracture  theory  to  materials  that  undergo 
some  plastic  deformation  during  crack  extension.  Later  in  1957,  Irwin  (5)  introduced  the 
concept  of  elastic  stress  intensity  factor  and  proposed  that  the  onset  of  crack  growth  in  a  cracked 
body  with  limited  plastic  deformation  at  the  crack  tip  is  when  the  elastic  stress  intensity  factor 
(K/)  reaches  a  critical  value  called  the  fracture  toughness.  Under  plane  strain  conditions,  the 
stress  intensity  factor  can  be  related  to  the  energy  release  rate  (G)  by  the  equation 


G  = 


(1) 


where  n  is  the  Poisson’s  ratio  and  E  is  the  elastic  or  Young’s  modulus.  There  are  numerous 
important  studies  on  dynamic  fracture  of  materials  that  were  executed  concurrently  with  the 
concepts  for  analyzing  cracked  solids  subjected  to  dynamic  loads.  For  more  information  on 
experimental  dynamic  fracture  research  of  engineering  materials,  the  reader  is  encouraged  to 
review  refs  6-15. 
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2.  Ductile  Fracture  of  Homogeneous  Materials 


In  general,  the  tenn  “ductile  fracture”  implies  that  during  a  failure  or  fracture  process,  the 
fractured  material  suffered  large  amounts  of  plastic  deformation.  This  is  somewhat  misleading 
since  cleavage  or  brittle  failure  sometimes  occurs  after  considerable  plastic  defonnation  in  some 
materials.  Usually  under  quasi-static  loading  conditions,  ductile  fracture  is  characterized  by  a 
cup-and-cone  failure  surface,  as  shown  in  figure  1 .  This  classic  form  of  failure  is  a  result  from  a 
complex  process  arising  from  the  stress  distribution  in  a  uniaxial  tensile  test  and  also  the  work¬ 
hardening  characteristics  of  the  material  being  tested.  There  are  other  types  of  failure  surfaces 
resulting  from  ductile  fracture,  such  as  the  double  cup  and  cone,  which  is  usually  observed  in 
ductile  materials  with  an  average  work  hardening  and  a  low  inclusion  density.  Another  failure 
surface  type  is  a  planar  surface  that  is  associated  with  high  strength  materials  with  high  inclusion 
content  and  a  low  work-hardening  rate. 


Figure  1.  A  classic  cup-and-cone 
failure  surface  resulting 
from  a  ductile  fracture  under 
quasi-static  load  conditions. 

How  these  macroscopic  failure  surfaces  formed  is  dependent  on  the  microscopic  fracture 
mechanisms  and  the  development  of  the  neck.  Usually,  a  neck  can  occur  when  the  increase  in 
tensile  stress  due  to  the  reduction  in  cross-sectional  area  as  a  result  of  plastic  defonnation  can  no 
longer  be  supported  by  the  concurrent  strain  hardening.  A  triaxial  state  of  stress  can  exist  in  the 
necked  region  after  a  neck  is  formed  due  to  the  nonunifonnity  of  the  necked  region.  This  triaxial 
stress  state  is  extremely  important  when  controlling  the  fracture  characteristics  of  a  material. 
Metallographic  studies  performed  on  the  necked  region  of  deformed  materials  have  shown  that 
during  the  early  stages  of  necking,  slip  is  the  dominating  factor.  However  for  large 
deformations,  voids  are  nucleated,  then  grow,  and  coalesce  to  form  cracks.  These  cracks 
ultimately  cause  material  failure.  The  nucleation,  growth,  and  coalescence  of  these  voids  are 
extremely  susceptible  to  stress  triaxiality. 


3 


2.1  Spallation  and  Fragmentation:  Void  Nucleation  and  Growth 

The  basic  definition  of  spall  fracture  is  when  two  plane  decompression  waves  subjected  to 
uniaxial  strain  conditions  collide  to  produce  a  region  of  tension  in  the  interior  of  a  material  body 
sufficient  to  cause  fracture.  Plate  impact  experiments  are  typically  used  for  studying  spall 
fracture.  As  shown  in  figure  2,  the  plate  impact  specimen  test  setup  mainly  consists  of  a  flyer 
and  a  target  plate.  The  flyer  plate  is  normally  launched  using  either  a  gun  system  with  the  aid  of 
a  sabot  or  an  explosive  toward  a  stationary  target  plate.  Upon  impact,  two  compression  waves 
are  generated  at  the  flyer-target  interface,  as  shown  in  figure  3,  and  they  travel  in  opposite 
directions.  Once  they  get  to  the  free  surface  of  both  the  flyer  and  target,  they  are  reflected  as 
decompression  waves  (or  rarefaction  waves).  Upon  collision  of  the  two  decompression  waves,  a 
region  of  high  tension  is  produced.  Strain  rates  of  104  to  106  s'1  are  typical  for  plate  impact  tests, 
and  the  resulting  stress  can  be  extremely  high  when  compared  to  the  quasi-static  strength  of  even 
the  strongest  materials.  Therefore,  the  material  fails  by  spallation.  It  is  worthy  to  note  that  the 
back  face  or  free  surface  velocity  of  the  target  is  monitored  using  a  Velocity  Interferometer 
System  for  Any  Reflector  (VISAR).  The  pull  back  velocity  obtained  from  the  velocity  profile  of 
the  back  face  velocity  is  used  for  estimating  the  spall  strength  of  the  material  in  question. 

Hopkinson  ( 16)  was  credited  with  the  earliest  observation  of  spall  fracture,  and  later  in  the  early 
1960s,  it  was  subsequently  studied  in  some  depth  by  Kolsky  (13)  and  then  Rinehart  and  Pearson 
(17).  During  the  1960s,  research  on  spall  fracture  became  more  intense  and  more  quantitative  as 
approaches  were  taken  which  led  to  theories  that  took  into  account  the  complexity  of  material 
behavior.  The  bulk  of  this  research  was  conducted  by  three  groups:  the  Stanford  Research 
Institute  (SRI),  Los  Alamos  National  Laboratory,  and  Sandia  National  Laboratory  (SNL).  The 
majority  of  this  section  will  focus  on  the  early  work  at  SNL  and  SRI.  Starting  around  the  mid 
1970s,  D.  P.  Dandekar  and  his  colleagues  at  U.S.  Army  research  establishments  all  around  the 
United  States  were  developing  novel  spall  testing  techniques  and  studying  the  response  of 
different  materials,  especially  ceramics,  to  support  novel  armor  designs.  Their  research  as  it 
relates  to  ballistics  is  reserved  for  further  review  by  the  authors. 


Figure  2.  Plate  impact  flyer-target  experimental  setup. 
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Figure  3.  A  simplified  distance-time  (Lagrangian)  plot  showing  the  propagation  of  waves  in 
the  flyer  and  target  after  impact. 


The  effects  of  the  amplitude  of  the  shock  compression  wave  that  precede  the  tensile  rarefaction 
wave  were  investigated  by  Stevens  and  Tuler  (18)  using  plate  impact  experiments.  They 
concluded  that  there  is  no  significant  effect  of  varying  amounts  on  the  dynamic  precompression 
on  the  spall  strength  of  1020  steel  and  6061-T6  aluminum.  This  simply  means  that  shock  wave 
strengthening  does  not  have  a  significant  effect  on  the  spall  strength  of  1020  steel  and  6061-T6 
aluminum.  However,  previous  experiments  done  by  Buchanan  and  James  (19)  showed  that  the 
spall  strength  of  mild  steel  increased  linearly  with  the  precompression  amplitude  for  a  pulse 
duration  that  was  not  constant.  Stevens  and  Tuler  (18),  on  the  contrary,  took  into  account  this 
notion  and  kept  the  pulse  duration  constant. 

After  reviewing  the  current  spall  criteria  at  the  time  and  systemizing  them,  Davison  and  Stevens 
(20)  introduced  the  concept  of  a  continuum  measure  of  spalling.  They  classified  the  spall  criteria 
into  classes,  “instantaneous  and  cumulative”  and  “local  and  non-local.”  If  the  spall  criteria  are 
dependent  only  on  the  current  values  of  the  field  variables,  then  the  class  is  “instantaneous”  and 
“cumulative”  if  the  spall  criteria  are  dependent  on  the  history  of  the  field  variables.  Furthermore, 
if  the  values  of  field  variables  at  the  spall  plane  are  used  in  the  detennination  of  damage,  then  the 
class  is  local  and  non-local  if  the  field  variables  at  distant  locations  have  a  role  on  the 
determination  of  damage.  The  theory  proposed  by  Davison  and  Stevens  (20)  was 
phenomenological  in  that  no  detailed  mechanisms  for  the  initiation  and  propagation  of  micro¬ 
failures  were  included.  However,  in  1973,  Davison  and  Stevens  (21)  proposed  a  more  detailed 
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theory  for  spall  related  damage  in  which  failure  occurs  by  the  initiation  and  propagation  of 
cracks.  This  was  applicable  to  alloys  exhibiting  brittle  behavior  subjected  to  dynamic  loads.  In 
1977,  a  comprehensive  report  was  published  by  Davison  et  al.  (22)  on  the  theory  of  ductile  spall 
damage  where  spalling  occurs  by  void  nucleation,  growth,  and  coalescence.  Figure  4  is  a 
schematic  representation  of  ductile  and  brittle  spall  fracture  resulting  from  a  projectile  impacting 
a  target. 


High  Strain  Rate  Deformation 


Dynamic  Fracture 


Sp  ailing  by  Brittle  Cracks 
Nucleation  Growth  and  Coalescence 


Spalling  by  Ductile  Voids 
Nucleation  Growth  and  Coalescence 


Adiabatic  Shear  Band 


Shock  Wave  Strengthening 


Figure  4.  Dynamic  deformation  system  (23). 

Over  the  past  two  decades,  several  theories  and  experimental  methods  have  been  developed  for 
predicting  the  behavior  and,  consequently,  the  fragmentation  resulting  from  explosively  driven 
cylinders  or  high-velocity  impacts.  Mott’s  theory  (24)  on  the  expansion  of  an  explosively  driven 
ring  and  its  resulting  fragmentation  was  expanded  by  Grady  (25).  Grady  calculated  the  fragment 
size  distribution  for  the  case  of  ductile  fracture  using  survival  statistics,  and  the  results  compared 
favorably  with  the  experimental  results  obtained  by  Wesenberg  and  Sagartz  ( 26)  for  an 
expanding  6061-T6  aluminum  cylinder.  A  general  expression  for  the  sizes  of  fragments  was 
cleverly  developed  by  Grady  based  on  an  energy  balance  approach.  The  basic  principles  of  the 
energy  balance  approach  are  that  the  interfacial  energy  generated  by  the  fragmentation  process  is 
balanced  by  the  local  inertial  or  kinetic  energy  of  the  material.  This  approach  was  first  applied  to 
the  dynamic  fragmentation  in  a  fluid,  as  shown  in  figure  5,  and  was  then  extended  to  brittle 
materials  by  incorporating  the  appropriate  fracture  mechanics  concepts.  The  instantaneous 
thennodynamic  and  kinematic  state  is  provided  by  the  density  p,  density  rate  p ,  and  temperature 
9,  which  may,  in  general,  be  functions  of  position.  Surface  tension  alone  resists  the  fracturing 
process.  It  is  intuitive  that  after  complete  fracture,  the  fragments  will  move  away  with  a  specific 
velocity  and  hence  kinetic  energy.  Because  the  fragments  continue  to  fly  away,  the  kinetic 
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energy  of  the  center  of  mass  of  each  fragment  is  maintained,  and  the  fragments  do  not  contribute 
to  the  generation  of  new  surfaces.  The  energy  terms  are  plotted  as  a  function  of  fracture  area  in 
figure  6,  and  the  fragment  size  is  determined  by  the  dE/dA  =  0,  which  is  the  minimum  of  the 
curve. 


Figure  5.  A  region  of  rapidly  expanding  fluid  that  will 
eventually  result  in  a  fragment  (25). 


Figure  6.  Kinetic  energy  vs.  surface  area  with  the 
minimum  as  the  fracture  area  (25). 
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At  equilibrium  (dE/dA  =  0),  the  fragment  surface  area  to  volumetric  ratio  is 


A  = 


"3p^V3 

spy 


(2) 


where  y  is  the  surface  energy  per  unit  area.  If  the  fragments  are  assumed  to  be  spherical  and  of 
equal  sizes,  then  the  fragment  diameter,  d,  is 

d  =  —  ■  (3) 

A 


Therefore, 


f 

d  =  6 

v 


spy 

3p2 


A 


J 


(4) 


The  elastic  stored  energy  has  to  be  considered  for  solids;  therefore,  the  equation  connecting 
fracture  mechanics  to  the  surface  tensions  is 


Y  = 


2pC2  ’ 


(5) 


for  which  KIC  is  the  fracture  toughness,  C  is  the  sound  speed,  and  p  is  the  density  of  the  material. 
If  equations  2,  3,  and  5  are  used  along  with  s  =  —  the  following  expression  for  the  fragment 

3p 


diameter  can  be  derived: 


d  = 


pCs 


2/3 


(6) 


Figure  7  shows  the  relationship  between  fragment  size  and  strain  rate  for  oil  shale.  There  is  good 
agreement  between  the  experimental  and  calculated  data. 


Grady’s  (25)  approach  could  be  generalized  by  including  the  crack-tip  plasticity  into  the 
calculations  as  explained  by  Hertzberg  (27).  The  effective  crack  length  is  computed  by 
including  some  fraction  of  the  plastic  defonnation  zone  ahead  of  the  crack,  and  for  plane  strain, 
the  plastic  zone  radius,  rv  is  estimated  to  be 


f  \2 

K 


6n 


vaw 


(7) 


where  K  is  the  stress  intensity  factor  and  ays  is  the  yield  stress.  Furthermore,  the  effective  stress 
intensity  factor  can  be  deduced  by  increasing  the  crack  length  ‘a’  by  ry  and  therefore 

Keff  =cr[7i(a  +  ry)}/2  .  (8) 
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Figure  7.  Comparison  between  calculated  and 
experimental  data  for  oil  shale  (25). 


By  substituting  equation  7  into  8,  the  effective  stress  intensity  factor  (Kefj)  can  be  determined  in 
terms  of  the  stress  intensity  factor  ( K ).  Because  the  onset  of  plasticity  is  determined  when  the 
stress  (cr)  is  equal  to  the  yield  stress  (crvv),  the  effective  stress  intensity  factor  (Keff)  can  be 
determined  to  be  1 .4 1  /f/c-  This  implies  that  equation  6  can  be  written  as 


d  = 


"dory 

pCs 


~ 6.3K1C  T 
pCs 


(9) 


which  is  the  diameter  of  the  fragment  in  terms  of  fracture  toughness,  density,  sound  speed,  and 
strain  rate. 

In  the  early  1970s,  the  SRI  played  a  vital  role  in  explaining  the  spall  failure  process  in  terms  of 
void  nucleation,  growth,  and  coalescence.  They  systematically  measured  cracks,  void  sizes,  and 
their  orientation  at  different  stages  of  the  spall  process  and  used  their  result  to  develop  the 
Nucleation  and  Growth  (NAG)  model.  They  developed  metallographic  techniques  that  are 
extremely  important  for  quantifying  parameters  describing  the  damage  accumulation  process. 

By  systematically  collecting  and  processing  data  for  a  number  of  materials,  Curran  et  al.  ( 28-36) 
laid  the  ground  work  for  the  fundamental  understanding  of  the  rate  of  nucleation  and  the  rate  of 
growth  of  micro-cracks  and  voids.  They  developed  the  equations  that  described  the  rate  of 
nucleation  N  and  the  rate  of  growth  R  (dR/dt)  to  be 

N  =  N()  exp[(cj  -  on0  )ct |  ]  (10) 

and 
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(11) 


R  = 


^  (j  —  a  ^ 

go 


4r) 


R 


where  N0  is  the  threshold  nucleation  rate,  ano  is  the  tensile  threshold  stress,  a,  is  the  stress 
sensitivity  for  nucleation,  <jg0  is  the  threshold  stress  for  growth,  and  r|  is  the  viscosity  of  the 
medium.  It  is  worth  noting  that  the  limiting  growth  rate  for  cracks  is  the  Rayleigh  wave  speed. 
From  equation  10,  it  is  clear  that  below  cno  ,  the  tensile  threshold  stress,  no  new  voids  are 
nucleated,  but  if  ano  is  surpassed,  the  nucleation  rate  increases  exponentially  with  the  tensile 
stress  a.  Equation  1 1  reveals  that  the  growth  rate  R  is  proportional  to  the  radius  of  the  void  or 
crack.  As  the  Rayleigh  wave  is  approached,  this  proportionality  breaks  down  for  a  crack.  It  has 
been  shown  that  for  spherical  voids  the  relative  void  volume  is  given  by 


y  =  .  8^AR03 


(a-(jgo)/4M 


exp 


3( 


a~Ggo\t 


4a 


-1 


(12) 


where  Ro  is  the  smallest  possible  void  radius.  This  was  assumed  to  be  1  pm  by  Barbee  et  al. 

(29).  A  one-  or  two-dimensional  (2-D)  Lagrangian  or  Eulerian  code  can  be  developed  to  predict 
relative  void  volumes  and  void  concentrations  as  a  function  of  distance  from  the  back  face  of  the 
target  specimen  once  the  parameters  of  equations  10,  11,  and  12  are  determined.  The 
NAG-FRAG  code  was  developed  by  SRI  using  the  approach  just  described. 

Interestingly,  the  mechanisms  for  ductile  failure  in  both  quasi-static  and  dynamic  load  conditions 
are  nucleation,  growth,  and  coalescence  of  micro-voids.  The  growth  of  voids  when  subjected  to 
quasi-static  loads  has  been  treated  theoretically  by  McClintock  (37,  38),  Rice  and  Tracey  (39), 
and  Needleman  (40).  McClintock  analyzed  cylindrical  voids  subjected  to  a  plane  strain 
condition  and  concluded  that  the  expansion  of  the  cylindrical  voids  eventually  leads  to  their 
coalescence.  His  analysis  also  shows  a  very  strong  inverse  dependence  of  fracture  strain  on  the 
tensile  stress  transverse  to  the  voids.  This  is  quite  important  to  the  spall  process,  where  the 
uniaxial  strain  condition  creates  strong  transverse  tensile  stresses  in  the  spall  zone.  And  as 
expected,  very  little  plastic  deformation  is  actually  observed,  hence  the  resulting  fracture  strain 
during  spallation  is  invariably  low.  Needleman  studied  the  interaction  between  the  stress  fields 
of  voids.  In  his  analysis,  double  periodic  arrays  of  parallel  cylindrical  voids  were  considered. 
Using  a  variational  method,  he  estimated  the  changes  in  void  shape  and  size  with  increasing 
overall  strain  and  was  also  able  to  track  the  development  of  the  plastic  zone  around  the  void. 

Using  a  first  approximation  of  the  variational  principle,  Rice  and  Tracey  (39)  analyzed  a 
spherical  void  of  radius  Ro,  as  shown  in  figure  8.  They  assumed  that  the  remote  strain  field 
consists  of  a  tensile  stretch  at  the  rate  s  with  respect  to  the  x3  direction  and  contractions  at  a  rate 

of  in  both  the  xi  and  X2  directions  (this  is  a  requirement  for  plastic  incompressibility).  Rice 

and  Tracey  found  the  velocity  field  to  be  in  the  fonn 
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Figure  8.  Spherical  void  in  a  remote  simple  tension 
strain  rate  field  (39). 


U:  =  £„XJC,  +  Du°  +  Euf  , 

1  lJ  J  1  1 


(13) 


where  D  and  E  are  constants,  utD  is  a  spherically  symmetric  volume  changing  field,  and  w/  is  a 

shape-changing  field  that  preserves  void  volume.  Because  of  incompressibility  and  spherical 
symmetry,  the  volume-changing  field  is  required  to  be 


•  D 
U;  =  8 


*0. 
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(14) 


Furthermore,  Rice  and  Tracey  (39)  considered  void  growth  in  a  general  remote  strain  rate  field 
with  high  stress  triaxiality  and  found  the  velocity  field  to  be 

ui  =  +  Du"  ,  (15) 


where 
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.  OO  •  00 
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Rice  and  Tracey  (39)  concluded  that  the  growth  rates  for  spherical  voids  are  significantly 
elevated  by  the  superposition  of  hydrostatic  tension  on  a  remotely  uniform  plastic  deformation 
field.  For  moderate  to  high  stress  triaxiality,  the  relative  void  growth  rates  are  amplified  over 
imposed  strain  rates  by  a  factor  depending  exponentially  on  the  mean  normal  stress. 


3.  Brittle  Fracture 


As  opposed  to  spallation/void  nucleation  and  coalescence,  the  fracture  of  materials  typically 
focuses  on  the  propagation  of  an  existing  crack  or  discontinuity.  The  highly  researched  field  of 
fracture  mechanics  spans  a  wide  range  of  rates  (to  include  fatigue,  quasi-static,  and  dynamic 
fracture),  types  (from  the  brittle  linear  elastic  fracture  mechanics  to  the  evolution  of  shear  cracks 
in  the  wake  of  shear  localizations),  and  materials  (including  glasses,  polymers,  metals,  adhesives, 
and  more).  In  this  section,  we  focus  on  high  strain-rate  brittle  fracture  of  materials. 

Linear  Elastic  Fracture  Mechanics  (LEFM)  is  the  classic  example  of  brittle  material  failure. 
Griffith’s  (2)  theory  of  fracture,  extended  to  elastodynamics,  revealed  some  interesting  results. 
Yoffe  (41)  considered  a  crack  of  constant  length  traveling  in  an  elastic  medium  at  a  fixed  speed; 
her  work  indicated  that  an  instability  in  the  problem  existed  as  the  speed  of  the  crack  approached 
the  Rayleigh  (surface)  wave  speed  (Cr).  Baker  (42)  extended  this  problem  to  a  transient  analysis 
and  demonstrated  a  similar  prediction  of  bifurcation  of  the  crack  when  its  velocity  approached 
the  Rayleigh  wave  speed.  Freund  (43)  offers  a  comprehensive  review  of  the  work  in 
elastodynamics  in  his  monograph,  to  include  the  analytical  solutions  for  modes  I,  II,  and  III 
cracks.  The  idea  of  a  “limiting  velocity”  sparked  many  experiments  and  assessments  of  crack 
velocities  with  respect  to  loading,  some  of  which  will  be  discussed  subsequently. 

Schardin  (6)  and  his  coworkers  utilized  a  multiple-spark  camera  in  an  attempt  to  observe  crack 
propagation  at  increasing  crack  velocities  in  various  types  of  glass.  They  confirmed  the 
existence  of  a  limiting  velocity  for  crack  propagation  and  showed  that  at  higher  driving  forces, 
the  crack  bifurcates  into  two  or  more  cracks  rather  than  increasing  in  speed.  However,  the 
limiting  velocity  that  Schardin  observed  was  substantially  lower  than  the  Rayleigh  wave  speed  in 
some  materials  it  was  as  low  as  30%  Cr.  Furthermore,  this  limiting  velocity  was  material 
dependent,  and  Schardin  proposed  that  this  was  an  additional  material  property  rather  than  an 
independent  mechanical  mechanism.  Kerkhof  (9)  utilized  a  separate  experimental  technique 
where  an  ultrasonic  transducer  input  a  small  amplitude  stress  wave  during  the  fracture  event, 
resulting  in  ripples  on  the  fracture  surface  whose  frequency  can  be  correlated  to  the  crack 
velocity.  Kerkhof  s  technique  confirmed  Schardin’s  findings,  and  by  examining  a  wider  variety 
of  glasses,  he  found  that  the  composition  of  the  glass  significantly  affected  the  limiting  crack 
velocity. 
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The  reduced  (less  than  the  Rayleigh  wave  speed)  limiting  velocity  has  been  the  subject  of  much 
discussion.  Ravi-Chandar  and  Knauss  ( 44-47)  conducted  an  extensive  study  of  dynamic  fracture 
in  a  brittle  polyester-based  polymer  (Homalite  100),  noting  that  the  surface  finish  of  the  fracture 
surfaces  shows  distinct  regimes  during  crack  propagation.  The  transition  from  “mirror”  to 
“mist”  to  “hackle”  showed  that  the  crack  demonstrated  an  increasing  tendency  toward 
branching/unsteady  propagation  as  the  crack  velocity  increased.  The  authors  proposed  that 
effects  at  the  micro-scale  played  a  significant  role  in  the  macro-scale  crack  propagation.  The 
formation  of  micro-cracks  ahead  of  the  macro-crack  tip  led  to  “discontinuous”  or  start/stop  crack 
motion  as  the  crack  jumped  between  micro-cracks,  and  these  micro-cracks  began  to  form  in 
more  substantial  numbers  as  the  crack  velocity  increased,  which  resulted  in  bifurcation  at  some 
critical  point.  The  formation  of  these  micro-cracks  is  inconsistent  with  the  linear  elastic  theory, 
explaining  the  ( 1 )  significantly  lower  limiting  velocity  than  predicted  by  traditional  theory  and 
(2)  difficulty  in  measuring  a  dynamic  fracture  toughness  (Kd)  that  predicts  the  behavior  of  this 
type  of  material.  In  contrast,  most  metals  undergo  yielding  near  the  crack  tip  without  the 
formation  of  the  micro-cracks.  The  linear  elastic  theory  can  be  modified  to  account  for  some 
localized  plasticity  (e.g.,  see  Freund  [43}),  and  the  resulting  theory  can  be  used  to  predict  a 
dynamic  fracture  toughness.  Experimental  efforts  have  shown  that  it  is  possible  to  measure  the 
dynamic  fracture  toughness  of  a  metal  as  a  function  of  crack  speed  (e.g.,  Pandolfi,  et  al.  [48}), 
and  that  theoretical/numerical  models  can  replicate  experimental  results. 

The  dynamic  propagation  of  mixed-mode  cracks  has  also  been  of  great  interest,  especially  given 
the  impact/penetration  scenario  as  a  driving  impetus.  Kalthoff  and  Winkler  (49)  introduced  the 
loading  cracks  by  edge  impact  (LECE1)  technique  to  examine  cracks  propagating  in  metals  for 
shear  (mode  II)  and  mixed  mode  (mode  I/II)  conditions.  Using  the  method  of  caustics,  Kalthoff 
and  Winkler  observed  that  failure  modes  can  transition  from  a  principal  stress  (tensile)  crack  to  a 
primarily  shear  (mode  II)  crack  by  increasing  the  driving  force  significantly,  implying  that  mode 
II  cracks  are  significantly  higher  energy  than  pure  mode  I  cracks  in  the  same  material. 
Furthermore,  for  some  materials  the  shear  crack  was  ductile  in  nature,  primarily  resulting  from 
the  formation  of  a  shear  localization  (adiabatic  shear  band)  and  the  resulting  void  nucleation  and 
coalescence.  Other  researchers  (e.g.,  Zhou  et  al.  [50])  have  observed  similar  behavior,  including 
the  formation  and  propagation  of  shear  bands,  in  metallic  materials.  Pure  mode  II  cracks  have 
exhibited  some  interesting  behavior  for  specific  systems.  Rosakis  et  al.  (51)  observed  mode  II 
cracks  propagating  along  weak  planes  in  a  brittle  polyester  resin  at  speeds  in  excess  of  the  shear 
wave  speed  of  the  bulk  material,  resulting  in  shear  shock  waves.  This  phenomenon  (detailed  in 
Rosakis  [52]),  particularly  of  interest  to  seismologists  and  earthquake  engineers,  allows  cracks  to 
propagate  at  extremely  high  speeds  in  high-energy  events  where  other  crack  modes  are 
prohibited  by  confinement. 
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Dynamic  brittle  fracture  has  a  rich  history  and  a  tremendous  amount  of  past  research — it  would 
be  impossible  to  reference  every  researcher  here.  Freund’s  (43)  monograph  and  the  book  by 
Ravi-Chandar  (53)  give  a  more  comprehensive  theoretical  discussion  of  dynamic  fracture 
research. 


4.  Numerical  Methods  in  Dynamic  Failure 


Accurate  simulation  of  the  dynamic  failure  of  materials  has  been  an  elusive  goal  in  the  material 
modeling  and  numerical  methods  communities  for  over  20  years.  Both  communities  have  a 
similar  problem — dynamic  failure  is  an  inherently  local,  discontinuous  change  in  the  behavior  of 
the  material  through  the  creation  of  new  surfaces.  For  material  modelers,  the  local  nature  of  the 
failure  makes  it  difficult  to  observe  in  the  laboratory,  and  the  length  and  timescales  of  dynamic 
failure  frequently  require  consideration  of  material  microstructure  and  defects.  The  discontinuity 
produced  by  failure  complicates  many  traditional  numerical  methods  (such  as  standard  finite 
element  fonnulations),  which  rely  on  continuous  state  variables  over  a  finite  discretization  for 
their  solutions.  Each  community  has  developed  several  interesting  methods  to  approach 
dynamic  failure. 

Computational  methods  for  simulating  dynamic  failure  can  usually  be  categorized  as  either 
continuum  or  discrete  methodologies.  Continuum  methods  are  typically  implemented  much  like 
constitutive  models  in  that  material  “failure”  (such  as  loss  of  deviatoric  or  total  strength)  is 
denoted  at  an  integration  point  of  a  finite  element  (or  other  numerical  method).  Discrete 
methods  make  an  attempt  to  denote  the  formation  of  new  surfaces,  either  through  mesh/particle 
re-mapping  or  through  the  introduction  of  additional  kinematics  that  accounts  for  the  formation 
of  the  surfaces.  Each  method  has  its  benefits  and  limitations,  sometimes  frequently  controlled  by 
implementation  as  much  as  physical  accuracy. 

Continuum  damage  mechanics  (CDM)  focuses  on  the  introduction  of  a  “damage”  variable  that 
describes  the  degradation  of  strength  and/or  modulus.  Introduced  in  1958  by  L.  M.  Kachanov 
(54),  CDM  has  spawned  a  tremendous  number  of  material/failure  models  ranging  from  plasticity 
to  brittle  failure  to  porosity  evolution.  CDM  assumes  that  the  irreversible  processes  occurring  on 
length  scales  below  the  one  of  interest  can  be  captured  using  a  set  of  state  variables  whose 
evolution  is  controlled  by  thennomechanical  phenomena  measurable  on  the  current  length  scale. 
Ductile  failure  models  such  as  the  Johnson-Cook  failure  (55),  Gurson  (56),  and  TEPLA  (57,  58) 
models  are  designed  with  CDM  in  mind;  the  length  scale  of  void  nucleation  and  growth  is  such 
that  explicit  representation  is  difficult  to  capture  with  realistic  computational  times.  Ballistic 
impact  on  brittle  materials  results  in  the  evolution  of  innumerable  micro-cracks;  in  an  attempt  to 
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model  this  behavior,  numerous  researchers  (e.g.,  Rajendran  and  Grove  [59],  and  Johnson  and 
Holmquist  [60])  have  used  CDM  methodologies  to  denote  the  loss  of  strength  as  the  micro-crack 
density  increases. 

CDM  is  also  used  extensively  in  the  composite  materials  community  because  of  the  multiple 
failure  modes  present  in  a  heterogeneous  material  with  a  small  internal  length  scale.  Most 
“multi-scale”  methods  utilize  a  fonn  of  CDM  to  transfer  information  from  the  meso-scale  to  the 
macro-scale;  for  example,  the  multi-continuum  theory  for  composite  materials  proposed  by 
Mayes  and  Hansen  (61)  transfers  damage  from  the  fiber  and  matrix  constituents  to  the  macro¬ 
scale  through  a  set  of  internal  state  variables.  One  principal  advantage  of  CDM  is  that  it 
accommodates  progressive  damage — certain  failure  modes/types  demonstrate  a  non-catastrophic 
loss  in  material  strength  (such  as  porosity  evolution),  which  is  difficult  to  capture  with  discrete 
methods  without  resorting  to  simulating  individual  voids.  CDM  is  also  straightforward  to 
implement — in  many  cases,  it’s  simply  adding  an  additional  subroutine  to  existing  material 
models.  However,  CDM’s  implementation  at  integration  points  also  leads  to  a  high  degree  of 
mesh  sensitivity  and/or  mesh  texturing  for  localized  failures  (such  as  brittle  cracking). 
Furthermore,  as  the  stiffness  of  the  element/particle/integration  unit  is  reduced  to  nearly  zero, 
distortions  can  become  unacceptably  large,  and  a  secondary  method  of  damage  (element 
deletion,  void  insertion,  etc.)  is  needed. 

Cohesive  zone  methods  have  gained  popularity  over  the  past  1 5  years  because  of  their 
straightforward  application  to  finite  element  analyses.  Cohesive  theories  of  fracture,  such  as 
those  proposed  by  Dugdale  (62)  and  Barenblatt  (63),  offer  a  different  theory  of  what  is 
happening  in  the  process  zone  near  a  crack  tip.  These  theories  assume  that  in  a  small  region 
ahead  of  the  crack  tip,  the  area  of  adhesion  decreases  progressively,  thus  resulting  in  a  lowered 
traction  along  the  “new”  crack  surface.  The  traction  vanishes  at  a  specified  opening 
displacement,  resulting  in  an  irreversible  energy  loss,  which  is  taken  to  represent  the  Griffith 
energy.  The  idea  of  a  process  zone  and  the  removal  of  the  singularity  present  in  linear  elastic 
fracture  mechanics  (and  thus  the  incorporation  of  an  inherent  length  scale)  make  cohesive  zone 
theory  ideally  suited  for  discretized  methods. 

Xu  and  Needleman  (64)  demonstrated  cohesive  zones  by  the  insertion  of  cohesive  surfaces 
between  all  of  the  elements  in  a  2-D  analysis,  assuming  an  effective  traction-displacement 
relation  between  the  cohesive  surfaces.  Transient  analyses  on  initially  cracked  plates  under 
impact  loading  showed  that  the  propagation  of  a  crack  and  its  ensuing  bifurcation  can  be 
qualitatively  captured  using  this  method.  Camacho  and  Ortiz  (65)  also  conducted  2-D  transient 
analyses  of  impact  but  utilized  the  adaptive  insertion  of  cohesive  zones — that  is,  a  traditional 
finite  element  analysis  with  full  connectivity  is  conducted  until  an  element  boundary  reaches  a 
specified  critical  traction,  at  which  time,  connectivity  is  modified,  a  cohesive  “element”  is 
inserted,  and  the  analysis  proceeds.  This  adaptive  insertion  of  cohesive  elements  prevents  some 
of  the  wave  dispersion  effects  seen  in  the  Xu  and  Needleman  approach  while  simultaneously 
reducing  the  computational  burden  by  considering  less  cohesive  zones.  The  method  has  since 
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been  extended  to  three  dimensions  (66).  The  cohesive  zone  methods  suffer  from  mesh 
dependency,  both  in  terms  of  element  size  and  path  dependency.  The  inherent  length  scale  of  the 
cohesive  theory  can  be  used  to  reduce  or  even  eliminate  the  mesh-size  dependency  by  computing 
the  process  zone  size — Falk  et  al.  (67)  consider  this  subject  in  depth — although  the  required 
element  size  may  produce  a  numerically  intractable  problem.  The  mesh  path  dependency, 
however,  cannot  be  fully  resolved  since  the  fracture  must  occur  along  element  boundaries, 
although  there  is  some  suggestion  that  random  meshing  using  triangles/tetrahedra  produces  a 
sufficient  number  of  crack  paths  to  reach  a  converged  solution  (68,  69).  The  cohesive  zone 
method  has  been  used  to  study  dynamic  fracture  along  weak  planes  (70),  penetration  into  brittle 
materials  (65),  fragmentation  of  brittle  materials  (71),  interface  crack  growth  (72),  and  various 
other  dynamic  failure  events. 

The  extended  Finite  Element  Method  (XFEM)  developed  by  Belytschko  and  his  coworkers  (see 
Moes  et  al.  [73])  has  also  proven  useful  in  the  simulation  of  discrete  failure  in  materials.  XFEM 
utilizes  enriched  nodes  (via  a  partition  of  unity  method)  to  allow  for  moving  discontinuities  in 
the  finite  element  mesh.  As  such,  crack/discontinuity  growth  is  not  limited  to  element 
boundaries,  and  the  solution  becomes  mesh-independent  as  the  overall  finite  element  solution 
converges.  The  nodal  enrichment  in  XFEM  is  problem/constitutive  dependent;  Moes  et  al.  (73) 
embedded  the  square -root  singularity  from  linear  elastic  fracture  mechanics  within  their 
formulation.  Obviously,  changes  in  constitutive  behavior  of  the  material  and/or  changes  in  the 
physics  of  failure  must  be  directly  incorporated  in  the  nodal  enrichment.  The  original  2-D  work 
was  extended  to  three  dimensions  in  Moes  et  al.  (74),  and  although  the  extension  is  nontrivial, 
the  authors  demonstrate  that  three-dimensional  (3-D)  XFEM  crack  modeling  is  possible  for 
elastic  and  inelastic  problems.  As  the  physics  for  the  cracking  phenomena  must  be  explicitly 
represented  in  the  XFEM  formulation,  XFEM  does  not  account  for  crack-branching  in  its 
standard  form;  however,  researchers  have  introduced  enrichments  that  allow  for  branched  and 
intersecting  cracks  (75). 

Song  et  al.  (69)  compare  XFEM  with  the  cohesive  zone  method  and  element  deletion,  noting  that 
both  the  XFEM  and  cohesive  zone  method  can  capture  dynamic  crack  propagation,  including 
bifurcation  in  mode  1  cracks.  Both  methods  require  small  elements  in  the  vicinity  of  the  crack  tip 
for  convergent  solutions — this  fact  is  dictated  by  the  small  size  of  the  physical  process  zone 
rather  than  any  numerical  deficiency.  Unlike  the  continuum  damage  techniques,  XFEM  and  the 
cohesive  zone  method  create  new  surfaces,  which  are  part  of  their  appeal,  but  they  also  create 
more  surfaces  for  the  numerical  contact  algorithms  as  well.  This  introduction  of  additional 
contact  surfaces  is  particularly  relevant  in  impact/penetration  problems  where  the  compressive 
state  of  stress  places  many  of  the  crack  faces  in  contact  with  one  another. 
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5.  Experimental  Techniques 


The  term  “dynamic  testing”  is  used  quite  loosely  to  cover  a  wide  range  of  strain  rates  and  testing 
equipment.  Although  one  would  wish  to  develop  material  properties  for  constitutive 
relationships  to  describe  a  material’s  behavior  across  a  wide  spectrum  of  strain  rates,  that  may 
not  usually  be  the  case  in  general.  This  is  to  be  expected  because  the  dominant  deformation 
mechanism  of  a  material  may  change  from  one  strain  rate  regime  to  another.  Figure  9  is  a  good 
summary  of  the  various  dynamic  characteristics  associated  with  mechanical  testing  (76). 
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Figure  9.  Dynamic  aspects  associated  with  mechanical  testing  (76). 

5.1  Plate  Impact 

Flere  we  discuss  gas  gun  driven  plate  impact  experiments  of  condensed  matter  (the  methodology 
described  here  is  generalized).  The  primary  objective  of  a  plate  impact  experiment  is  to  generate 
a  planar  shock  compression  wave  in  a  uniaxial  state  of  strain  in  the  target  material.  This  is 
normally  attained  by  the  use  of  gas  guns  (single  stage,  two  stage,  powder,  or  electromagnetic), 
explosives,  or  lasers,  depending  on  the  strain  rate  of  interest.  With  the  use  of  a  gas  gun,  a  flyer 
plate  mounted  on  a  sabot  is  driven  into  a  target,  resulting  in  a  symmetric  or  non-symmetric 
impact,  depending  on  the  flyer  and  target  material.  If  the  flyer  and  target  material  are  identical, 
then  the  impact  is  said  to  be  symmetric;  otherwise,  it  is  non-symmetric.  A  series  of  charged  pins 
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is  used  to  determine  the  flyer  velocity  prior  to  impacting  the  target.  For  a  test  to  be  considered 
acceptable,  the  impact  between  the  flyer  and  the  target  must  be  parallel  to  within  2  mrad.  The 
lateral  dimensions  of  the  specimens  are  chosen  so  that  uniaxial  strain  condition  is  maintained 
throughout  the  time  of  interest,  during  which  free  surface  velocity  measurements  are  made.  In 
other  words,  free  surface  velocity  measurements  are  made  prior  to  the  arrival  of  the  radial  release 
waves.  Upon  impact,  a  two-structure  shock  compression  wave  is  generated  at  the  impact 
interface,  traveling  in  opposite  direction  as  shown  in  figure  10  (a),  the  time-distance  or 
Lagrangian  plot.  The  two-structure  shock  compression  wave  comprises  an  elastic  wave  and  a 
slower  plastic  wave.  At  time  ti,  the  first  time-of-arrival  signal  for  the  elastic  wave  is  detected  by 
the  Velocity  Interferometer  System  for  Any  Reflector  (VISAR  [77])  at  the  free  surface  of  the 
target  and  shortly  thereafter  at  time  t2  that  of  the  plastic  wave  is  detected. 


Figure  10.  (a)  Lagrangian  or  time-distance  diagram  and  (b)  the  velocity-time 
history  of  a  two-structure  wave. 

Both  waves  are  reflected  as  rarefaction  waves.  Similarly,  the  two-structure  shock  compression 
wave  in  the  flyer  specimen  is  reflected  at  the  free  surface  as  a  rarefaction  wave.  Assuming  a 
symmetric  impact,  both  elastic  and  plastic  rarefaction  waves  from  the  flyer  free  surface  are 
completely  transmitted  through  the  impact  interface  into  the  target  material.  When  the 
rarefaction  waves  from  both  free  surfaces  interact,  a  region  of  almost  pure  hydrostatic  tension  is 
generated,  and  stresses  that  are  greater  than  the  threshold  stress  required  for  damage  initiation  are 
produced,  as  shown  in  figure  10  (a).  Within  this  region,  voids  nucleate,  grow,  and  coalesce  to 
form  well-defined  spall  damage  in  ductile  materials  and  can  result  in  material  failure.  A  good 
review  of  the  plate  impact  technique  has  been  compiled  by  Fowles  et  al.  (75). 

The  discovery  of  lasers  has  led  to  the  invention  of  optical  interferometery,  which  currently  is  the 
most  powerful  tool  in  the  determination  of  the  surface  motion  of  shocked  condensed  matter.  A 
schematic  representation  of  the  velocity-time  history  of  a  plate  impact  experiment  is  shown  in 
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figure  10  (b).  Acquiring  the  velocity-time  history  is  of  the  utmost  importance  in  spall  fracture 
experiments.  Useful  experimental  information  such  as  the  Hugoniot  Elastic  Limit  (HEL),  plastic 
rise  time,  peak  pressure  (related  to  the  maximum  velocity,  Vmax),  pull-back  velocity  (AUfs),  etc., 
can  be  extracted  from  the  velocity-time  history.  The  HEL  is  defined  as  the  point  at  which  a 
solid,  loaded  in  compression  under  the  constraint  of  uniaxial  strain,  can  no  longer  support  elastic 
distortion  and  begins  to  flow  through  plastic  or  cataclastic  (crushing)  fracture  processes  (79). 
Also,  the  spall  strength  is  directly  related  to  the  pull-back  velocity  by  the  relationship 

*  1 

°  =2poco  AUjs,  (17) 

where  a*  is  the  spall  strength,  p0  is  the  density  of  the  material,  and  C0  is  the  bulk  sound  speed. 
Acquiring  the  velocity-time  history  from  a  plate  impact  experiment  requires  a  VISAR.  The 
VISAR  has  been  one  of  the  most  important  breakthroughs  in  experimental  shock  physics  of 
condensed  matter.  The  VISAR  is  an  extension  of  the  Michelson  Interferometer,  and  the  working 
principles  are  based  on  Doppler  shift,  i.e.,  due  to  the  movement  of  the  free  surface  of  the  target 
specimen,  the  system  produces  interference  fringe  shifts  that  are  proportional  to  the  Doppler  shift 
of  the  laser  light.  Because  of  the  importance  of  VISAR  in  experimental  shock  physics,  its 
operation  is  discussed  in  the  following  paragraphs.  A  schematic  of  the  VISAR  system  is  shown 
in  figure  1 1 .  The  impact  between  the  flyer  and  target  takes  place  inside  a  vacuum  chamber,  as 
shown  in  the  figure.  A  laser  beam  is  focused  on  the  back  free  surface  of  the  target,  and  the 
reflected  light  is  collimated  and  then  split  into  two  beams  by  a  10-90  beam  splitter,  i.e.,  10% 
reflected  and  90%  transmitted.  The  10%  reflected  beam  is  directed  to  a  Beam  Intensity  Monitor 
(BIM)  for  monitoring  and  diagnostic  purposes.  The  90%  transmitted  beam  is  further  split  into  a 
reference  leg  and  a  delay  leg  by  a  50-50  beam  splitter.  The  delay  leg  is  achieved  by  transmitting 
the  beam  through  a  glass  etalon  (fused  silica)  and  then  shifted  by  a  X/8  wave  plate.  The  purpose 
of  the  delay  is  to  achieve  a  90°  phase  difference  between  the  two  beams.  The  physical  distances 
traveled  by  the  reference  and  delay  legs  are  made  to  be  equivalent  to  within  a  few  thousandths  of 
an  inch  in  order  to  obtain  quality  fringe  patterns.  Both  legs  are  then  separated  into  S  and  P 
components  of  light  by  a  Polarizing  Beam  Splitter  (PBS)  and  then  recombined  to  form  two  90° 
out-of-phase  signals  that  are  sinusoidal  in  nature.  These  fringes  are  shifted  whenever  the  free 
surface  of  the  target  is  accelerated,  and  surface  velocity  V(t )  is  related  to  the  total  fringe  count 
F(t)  by  the  equation 

F-(,-It)  =  V - - •  (18) 

2  2t(1  +  6)(1  +  Hv/v) 

where  X  is  the  wavelength  of  the  laser  beam,  5  is  a  correction  factor  due  to  the  wavelength 
dependence  of  the  refractive  index  on  the  etalon,  and  Av/v  is  an  optical  correction  term. 
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Figure  1 1.  A  schematic  representation  of  the  VISAR  system. 

5.2  Dynamic  Fracture 

Experimental  techniques  to  elucidate  dynamic  fracture  have  taken  many  forms  over  the  past  50 
years,  but  all  focus  on  the  ability  to  discern  the  crack  speed  and,  more  recently,  the  state  of 
stress/strain  in  the  vicinity  of  the  crack  tip.  Techniques  can  generally  be  classified  as  optical  and 
non-optical,  with  current  research  trending  toward  optical  methods  (due  to  their  ability  to  capture 
full-field  information). 

Non-optical  methods  for  dynamic  fracture  can  be  further  subdivided  into  stress  wave/crack 
interaction  and  electrical  resistance  methods.  Wallner  (80)  discovered  that  during  crack 
propagation,  the  interaction  between  a  propagating  crack  and  shear  waves  produced  during  the 
fracture  process  produced  distinct  ridges  in  the  fracture  surface.  In  a  postmortem  analysis,  one 
can  use  these  ridges  (along  with  the  shear  wave  speed  in  the  material)  to  estimate  the  crack 
velocity. 

Ravi-Chandar  (53)  gives  an  excellent  review  of  the  technical  details  of  these  “Wallner  lines”  and 
how  they  can  be  used  to  measure  the  crack  speed.  Kerkhof  (9)  extended  this  method  by 
introducing  an  ultrasonic  transducer  as  a  source  of  shear  waves;  now  the  creation  of  the  Wallner 
lines  is  governed  by  a  user-controlled  frequency,  thus  appearing  as  smooth  “ripples”  on  the  crack 
surface  with  a  specified  temporal  frequency.  Their  spatial  frequency  can  then  be  used  to 
determine  the  crack  speed  with  a  high  degree  of  accuracy.  This  so-called  “stress  wave 
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fractography”  is  discussed  in  detail  in  Richter  and  Kerkhof  (81).  Kerkhof  (9)  used  this  method 
to  find  the  speed  of  cracks  in  glasses  with  a  1-MHz  transducer,  allowing  the  accurate 
determination  of  crack  speeds  up  to  2000  m/s. 

Electrical  resistance  methods  rely  on  sets  or  grids  of  fine  electrical  wires  applied  to  the  material 
of  interest.  As  the  crack  passes,  each  wire  is  broken,  providing  a  direct  crack  position  as  a 
function  of  time.  Many  researchers  have  utilized  this  technique — possibly  the  first  being 
Dulaney  and  Brace  (82) — and  the  grids  are  now  commercially  available.  Thin  films  coupled 
with  alternating  currents  can  also  be  utilized  to  determine  crack  velocities  (83)  at  a  high 
sampling  rate  when  great  care  is  taken  to  prepare  the  thin  film  and  measure  its  resistance  changes 
as  the  crack  propagates. 

Optical  techniques  utilizing  high-speed  cameras  are  some  of  the  most  prominent  dynamic 
fracture  instrumentation  techniques.  Crack  propagation  speeds  in  excess  of  1  km/s  require  the 
use  of  cameras  with  very  high  framing  rates  to  capture  fracture  processes  in  situ.  Schardin  (6) 
developed  a  spark  camera  (now  frequently  referred  to  as  the  Cranz-Schardin  camera)  that  utilizes 
a  set  of  light  sources  and  lenses  that  produces  separate  images  for  the  pulse  of  each  light  source; 
thus  framing  rates  approaching  1  million  frames  per  second  can  be  obtained,  as  the  exposure 
time  is  independent  of  the  framing  rate.  Lu  and  Liu  (84)  offer  a  review  of  the  construction  of 
Cranz-Schardin  cameras,  discussing  the  original  spark-based  system  as  well  as  a  modem  Light 
Emitting  Diode  (LED)  system.  Rotating  mirror  cameras  can  improve  the  temporal  frequency, 
approaching  5  million  frames  per  second,  while  still  retaining  excellent  image  quality.  Gated 
intensified  cameras  utilize  beam  splitting  optics  and  Charged  Coupled  Device  (CCD)  arrays  to 
generate  extremely  high  framing  rates  (>200  million  frames  per  second)  for  a  limited  number  of 
frames  at  very  short  exposure  times  (<5  ns)  with  reduced  image  quality  as  compared  to  rotating 
mirror  cameras.  Recently,  several  manufacturers  have  developed  single  CCD  video  cameras 
capable  of  framing  rates  in  excess  of  1  million  frames  per  second  with  excellent  image  quality 
for  dozens,  if  not  hundreds,  of  frames.  Regardless  of  the  camera  type,  several  imaging 
techniques  have  been  utilized  by  researchers  to  gain  insight  into  dynamic  fracture. 

Direct  high  speed  photography  grants  some  information  about  the  propagation  of  cracks  in 
transparent  materials,  but  other  techniques  are  often  employed  that  give  additional  information. 
The  shadowgraph  method  frequently  referred  to  as  the  method  of  caustics  (85)  utilizes  the  local 
slope  change  (opaque,  reflective  materials)  or  the  thickness  change  (transparent  materials)  near 
the  crack  tip  to  produce  a  shadow  “spot”  that  precisely  locates  the  crack  tip  during  a  dynamic 
event.  Detailed  analyses  by  Kalthoff  (86)  show  that  one  can  use  the  size  and  shape  of  this  spot  to 
determine  the  relative  magnitude  of  the  loading  and  the  mode  mixity  of  the  crack,  respectively. 
An  example  of  an  LECEI  test  for  PMMA  is  shown  in  figure  12  (87). 
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Figure  12.  Method  of  caustics  applied  to  edge-impacted  pre-cracked  PMMA  specimen,  resulting  in 
mode-II  dominant  cracks  (87). 


Tippur  et  al.  (55)  introduced  an  interferometric  technique  that  measures  gradients  of  the  out-of- 
plane  surface  displacements,  naming  the  technique  coherent  gradient  sensing  (CGS).  This 
method  has  proven  to  be  quite  useful  for  dynamic  crack  propagation;  it  is  relatively  insensitive  to 
vibrations,  entirely  insensitive  to  rigid  body  motions,  and  the  resulting  fringe  patterns  can 
distinguish  the  mode  mixity  of  the  crack  through  comparison  to  analytical  solutions.  The  full- 
field  information  allows  critical  examination  of  crack  bifurcation,  including  shifts  in  mode- 
mixity  or  failure  type  (such  as  a  shift  to  adiabatic  shear  banding  [59]).  Furthermore,  a  slight 
change  in  experimental  setup  allows  the  technique  to  be  applied  to  either  transparent  or 
specularly  reflective  materials.  Figure  13  shows  the  optical  setup  for  reflective  CGS;  details  of 
the  analysis  technique  can  be  found  in  Tippur  (90).  The  technique  uses  a  coherent  light  source 
and  a  set  of  carefully  aligned  Ronchi  rulings  to  create  a  series  of  fringe-based  images;  a  lens  is 
used  to  focus  each  image  to  a  small  enough  size  such  that  an  aperture  can  isolate  the  image  of 
interest.  The  +1/-1  images  (the  first  image  above  and  below  the  optical  axis)  produce  the 
gradient  of  the  out-of-plane  displacement  in  the  direction  of  the  Ronchi  rulings;  obviously,  the 
Ronchi  rulings  can  be  rotated  to  give  different  gradients  as  desired. 
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Figure  13.  Diagram  of  coherent  gradient  sensing  experimental 

setup  for  reflective  materials  (from  Rosakis  et  al.  [97]). 


CGS  has  been  successfully  employed  to  study  the  fracture  behavior  of  metals  (48),  metallic 
glasses  (92),  fiber-reinforced  composites  (93),  and  polymers  (94);  indeed  the  method  is 
applicable  to  any  material  with  sufficient  local  deformation  at  the  crack  tip.  CGS  is  difficult  to 
employ  for  materials  with  extremely  small  strains  at  fracture  (e.g.,  borosilicate  glass) — the  small 
strain  gradients  produce  too  few  fringes  to  reliably  measure.  While  not  a  limitation  of  the 
technique,  the  amount  of  light  required  for  high-speed  photography  dictates  a  significantly 
powerful  laser  with  a  long  coherence  length  (argon-ion  and  diode-pumped  solid  state  lasers  of 
approximately  10  W  have  been  used  with  exposure  times  of  approximately  0.5  ps).  Other  optical 
techniques  have  been  utilized  to  measure  deformations  of  the  sample  in  an  attempt  to  monitor 
dynamic  crack  propagation.  Digital  image  correlation  (DIC)  (95)  utilizes  a  random  speckle 
pattern  and  one  or  more  high-speed  cameras  to  track  the  motion  of  the  speckles  and  produce  a 
deformation  field  of  the  sample.  A  single  camera  can  track  in-plane  deformations;  adding  a 
second  with  an  appropriate  angle  can  possibly  allow  3-D  tracking.  DIC  has  been  very 
successfully  employed  in  tracking  strain  behavior  of  materials  and  defonnations  of  structures; 
crack  propagation  tests  the  method  because  of  the  very  small  opening  displacements,  particularly 
near  the  crack  tip.  Direct  interferometric  measurement  of  the  out-of-plane  displacements  has 
been  employed  by  numerous  researchers  (e.g.,  Kalthoff  [57])  for  dynamic  crack  propagation. 
Given  that  this  measurement  relies  on  the  wavelength  of  the  coherent  light  source  (frequently 
approximately  500  mn),  the  resulting  fringe  pattern  captures  very  small  deformations.  While 
quite  useful,  setup  and  post-processing  can  be  difficult  as  the  method  is  sensitive  to  initial 
alignment  and  rigid  body  rotations. 
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6.  Concluding  Remarks 


The  understanding,  characterization,  and  simulation  of  dynamic  failure  of  materials  are  critical 
for  military  weapons  research;  areas  such  as  warhead  fragmentation,  penetration  of  brittle 
armors,  and  projectile  breakup  are  controlled  by  material  failure  and  its  kinetics.  This 
knowledge  also  impacts  civilian  applications,  such  as  fragmentation  of  rocks  for  mineral  and  oil 
extraction,  protection  of  space  vehicles  from  cosmic  debris  impact,  barricade  design  for  civil 
structures.  The  presented  review  shows  the  advancements  in  dynamic  failure  research  including 
failure  modes  ranging  from  spallation  to  brittle  fracture  over  the  past  50  years.  Experimental 
techniques,  theoretical  models,  and  computational  methods  for  dynamic  failure  were  discussed, 
giving  a  historical  state-of-the-art  perspective  on  work  performed  in  this  area. 

While  there  is  an  abundance  of  work  in  the  field,  several  key  topics  resonated  with  the  authors 
for  the  focus  of  future  work.  Experimental  methods  are  trending  toward  higher  resolution  in 
both  space  and  time;  full-field  and/or  penetrative  (X-ray/proton  radiography)  imaging  in  situ  can 
give  some  indication  of  the  kinetics  of  dynamic  failure  that  is  not  apparent  from  postmortem 
examination.  These  experiments  can  provide  some  needed  verification/validation  for  both  old 
and  new  models  for  dynamic  failure.  Numerical  modeling  techniques  for  dynamic  failure  are 
becoming  more  sophisticated  (e.g.,  extended  finite  element  and  cohesive  zone  methods)  rather 
than  the  traditional  continuum  damage  mechanics  approaches.  This  trend  should  help  material 
modelers  better  represent  failure  and  failed  material,  as  there  is  now  a  numerical  representation 
of  the  evolving  surfaces — a  critical  step  to  correctly  modeling  wave  propagation  and  local  stress 
states,  ft  is  currently  unknown  how  well  the  current  trend  of  multi-scale  modeling  will  impact 
dynamic  failure;  however,  it  could  potentially  help  link  smaller-  scale  processes,  such  as 
dislocation  dynamics,  grain-to-grain  interface  failure,  and  pre-existing  microstructural  flaws  to 
the  nucleation  and  propagation  of  macro-scale  cracks.  Invariably,  the  materials  science  and 
mechanics  communities  will  need  close  collaboration  to  reach  a  fundamental  understanding  of 
the  dynamic  failure  of  materials,  a  difficult  but  rewarding  goal  to  achieve. 
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